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Abstract 

The stationary, axisymmetric reduction of the vacuum Einstein equations, the so- 
called Ernst equation, is an integrable nonlinear PDE in two dimensions. There 
now exists a general method for analyzing boundary value problems for integrable 
PDEs, and this method consists of two steps: (a) Construct an integral represen- 
tation of the solution characterized via a matrix Ricmann-Hilbcrt (RH) problem 
formulated in the complex fc-plane, where k denotes the spectral parameter of 
the associated Lax pair. This representation involves, in general, some unknown 
boundary values, thus the solution formula is not yet effective, (b) Characterize 
the unknown boundary values by analyzing a certain equation called the global 
relation. This analysis involves, in general, the solution of a nonlinear problem; 
however, for certain boundary value problems called linearizable, it is possible to 
determine the unknown boundary values using only linear operations. Here, we 
employ the above methodology for the investigation of certain boundary value 
problems for the elliptic version of the Ernst equation. For this problem, the 
main novelty is the occurence of the spectral parameter in the form of a square 
root and this necessitates the introduction of a two-sheeted Ricmann surface for 
the formulation of the relevant RH problem. As a concrete application of the 
general formalism, it is shown that the particular boundary value problem cor- 
responding to the physically significant case of a rotating disk is a linearizable 
boundary value problem. In this way the remarkable results of Neugebauer and 
Meinel are recovered. 
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1 Introduction 

The stationary axisymmetric vacuum Einstein equations can be reduced to a single 
nonlinear PDE, the so-called Ernst equation. The elliptic version of this equation is 
given by 

^ (f PP + /« + - p f P \ = f P + fl P > 0, C G M, (1.1) 

where /(/?, £) is a complex-valued function, called the Ernst potential and bar denotes 
complex conjugation. 

A unified method for analyzing boundary-value problems (BVPs) for linear and 
integrable nonlinear PDEs in two dimensions was introduced in [6] and developed 
further by several authors |Il|2|iaH|[aiSllIDliniII21IiaiIlinil221|231. This method 
consists of two novel steps: (a) Construct an integral representation of the solution by 
performing the simultaneous spectral analysis of both parts of the Lax pair (this is to 
be contrasted with the inverse scattering transform method where one only performs 
the spectral analysis of the t- independent part of the Lax pair), (b) Characterize 
the unknown boundary values by analyzing the so-called global relation. Step (a) 
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characterizes the solution in terms of a Riemann-Hilbert (RH) problem involving all 
boundary values. Thus, in order to obtain an effective solution, it is necessary to 
implement step (b) . The complexity of this step depends on the particular boundary 
value problem under consideration. For example, for linear evolution PDEs and for 
some simple BVPs for linear elliptic PDEs, the unknown boundary values can be 
obtained in closed form. On the other hand, for nonlinear evolution PDEs on the 
half-line, the unknown boundary values, in general, are characterized through the 
solution of a nonlinear Volterra integral equation; however, for particular boundary 
conditions called linearizable, step (b) can be solved in closed form. 

Here, we implement this methodology to the elliptic version of the Ernst equation 



(1.1). For this problem, the main novelty is the occurence in the Lax pair of the 
spectral parameter in the form of a square root, which necessitates the introduction 
of a two-sheeted Riemann surface. A most interesting feature of the Ernst equation 
is the existence of a large class of linearizable BVPs. Let T> denote the exterior of a 
finite disk of radius po > (see Figure |L~T] ) , i.e. 

V = {(p,() € R 2 \p > 0}\{(p,0) G M 2 |0 < p < po}, 
and let f'(p,C) denote the Ernst potential in coordinates corotating with angular 



velocity O (see subsection 5.2). Then the following BVPs are linearizable: 



(A) (/' = on a finite disk) 



(i) f satisfies (1.1) in V. 



(ii) f(p,C) = f(Pi — C) (equatorial symmetry). 

(Hi) f(p, C) — > 1 as p 2 + ( 2 — > oo (asymptotic flatness). 

(iv) f p (+0, C) = for all £ ^= (regularity on rotation axis). 

(v) f' p (p, ±0) = for < p < po (constant Dirichlet boundary conditions on 
the disk). 

(B) (f c = on a finite disk) 

The solution / satisfies (i)-(iv) of (A) but (v) is replaced with 

(v 1 ) f'^(p, ±0) = for < p < po (vanishing Neumann boundary conditions on 
the disk). 

(C) (f' p = or = on an infinite disk) 

The solution / satisfies the same BVPs as in (A) and (B), respectively, but the 
domain D is replaced with the exterior V of an infinite disk stretching from 
po > to infinity, i.e. 

V = {(p,()\p>0}\{(p,0)\p o <p}. 

In the case of /' = 0, the condition (Hi) of asymptotic flatness must be replaced 
with a condition compatible with the constant value of /' on the disk. 
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Figure 1.1 The exterior domain V of a finite disk of radius pq. 



The physically significant problem (A) has been studied extensively in the pioneer- 
ing work of Neugebauer and Meinel [18 -[20 (see also [21J). It is remarkable that these 
authors were able to solve this problem without the guidance of any general method. 
We can now revisit this problem with the advantage of having at our disposal the 
general methodology of [6] and this, we hope, makes the relevant construction easier 
to motivate and simpler. 

It was emphasized in [6] that the new methodology provides a new approach to 
solving linear BVPs. Actually, if a given linear BVP can be solved by the new method, 
then the corresponding nonlinear problem can be solved following conceptually similar 
(but analytically more complicated) steps. For this reason we will first investigate 
the linearized version of equation (|1.1[) . 



1.1 Organization of the paper 

In section [2] we analyze the axisymmetric Laplace equation, or equivalently the static 
version of the Ernst equation. In section [3] we implement step (a) of the new method 
for the exterior finite disk domain depicted in Figure [i~T} This yields an expression for 
the Ernst potential / in terms of the solution of a RH problem which involves certain 
spectral functions defined in terms of both the Dirichlet and Neumann boundary 
values on the disk, see Proposition 3.4 In section [4] we analyze the consequences of / 



being equatorially symmetric and derive the global relation. In section [5] we analyze 
the global relation for the particular case of the BVP specified in (A). The BVPs 
formulated in (B) and (C) will be analyzed elsewhere. 
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2 The axisymmetric Laplace equation 



The linearized version of equation is the axisymmetric Laplace equation given 
by 

f PP + fa + -f P = o, p > o, C g k. 



Instead of considering this equation as the linearized approximation of (1.1), it is 



convenient to view it as the exact formulation of ( 1.1 ) in the case of a static spacetime. 



A static (as opposed to stationary) space time corresponds to a real-valued Ernst 
potential / = e 2U , C/£R, and in this case (1.1) reduces to the following equation for 
the real- valued function U (p, £) : 



u PP + u a + -u p = o, 



1 

—I 

p 



p > 0, (Gi. 

The boundary conditions on the rotation axis and at infinity are 

C/ p (+0,C) = for 

U(p, C)->0 as /3 2 + C 2 ^oo. 



(2.1a) 



(2.1b) 
(2.1c) 



2.1 Lax pair 

It is convenient to work with a complex variable z = p+iC and write U(z) for ?7(p, C)- 
For convenience of notation we will suppress the dependence on z , i.e. in general f{z) 
will denote a function depending on both z and z. 



Equation (2.1a) admits the following Lax pair formulation for the scalar function 
4>(z,k): 

U z (z,k) = XU z (z), 
U- z {z : k) = \U- z (z), 



where A = \(z, k) is defined by 



k — iz\ 



k + iz 



and k E C is a spectral parameter. We write (2.2) in the differential form 



W, 



where W = W(z, k) is the one-form 



1 



W =\U z dz + -U s dz 



X 

l + \)u p + i(±-\)u c 



dp + 



+ X)U, 



(2.2) 



A U D 



(2.3) 
(2.4) 
(2.5) 
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2.2 A two-sheeted Riemann surface 



According to the methodology introduced in [7j, the solution of the so-called direct 



problem involves the construction of a solution <p(z,k) of equation (2.4), which is 



bounded for all k G C. For a polygonal domain, this can be achieved by integrating 



(2.4) along a contour starting at a corner of the domain. In particular, the corner ioo 
yields the solution 

<p(z,k) = [ W(z',k). (2.6) 



However, the definition (2.3) of A involves a square root, hence the value of the right- 



hand side of (2.6) will depend on the path of integration since the path affects the 



choice of the branch of the relevant square root. In order for the right-hand side 



of (2.6) to be independent of the path of integration (so that (f> is well-defined), we 
introduce the following genus Riemann surface S z \ For each value of z — p -\- S z 
consists of the set of points (A, k) G C 2 such that 

, o k %z 

A = -, — -. 

k + iz 

We introduce a branch cut in the complex /c-plane from — iz to iz and, for k G C, we 
let k + and k~ denote the corresponding points on the upper and lower sheet of S z , 
respectively. By definition, the upper (lower) sheet is characterized by A — » 1 (A —* 
— 1) as k — > oo. We compactify S z by adding the two points oo + and oo~. For each 



z, 4>(z, •) is a map S z — > C. The definition (2.6) of (f>(z, k) is made precise by choosing 



k to lie on the lower sheet initially when z' = too. This picks one of the two branches 
of the square root defining X(z', k). We then use analytic continuation to follow this 
branch throughout the integration. This defines the integrand unambiguously. At 
the end of the integration, i.e. when z' approaches z, k will lie either on the upper 
or on the lower sheet of S z ; we denote the corresponding two values of 4> by (f)(z, k + ) 
and (p(z,k~), respectively. 

2.3 The direct problem 

In order to define a function (f>(z, k) for all k on S z , starting from z' = ioo, we integrate 



with respect to the contours 71 and 72, see graphs on the left of Figure 2.1 It turns 
out that for £ > 0, the function <p defined with respect to 71 and 72 lives on the lower 
and upper sheet of S z , respectively, i.e. 

<t>(z,k~)= ! W(z',k), 4>(z,k + )= ! W(z',k). (2.7) 

J 7! J 72 

Indeed, as z' moves along 71 and 72, the endpoints of the branch cut in the fc-plane 



move along the dotted curves of Figure 2.1. Note that for z = iC, on the C-axis, the 



Riemann surface S z degenerates and consists of two disjoint copies of the complex 
k- plane; we have A = 1 for all k + on the upper sheet and A = — 1 for all k~ on the 
lower sheet. Therefore, when integrating along the £-axis, the branch cut disappears. 

For the movement of the branch cut we note that the branch points for z' = p'+iC' 
occur at (£' — ip' , + ip'). Any fixed k on the lower sheet will remain on the lower 
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Figure 2.1 The graphs on the left show the integration contours 71 and 72 used to 
define (f>(z, k~) and 4>(z, k + ) for ( = Im z > 0. The graphs on the right illustrate how 
the endpoints of the branch cut in the complex k-plane move as these contours are 
traversed. For z' on the (-axis the Riemann surface degenerates and the branch cut 
disappears. 
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sheet throughout an integration along 71, because k never crosses the branch cutQ 
On the other hand, as the integration along the large semicircle of 72 is performed, 
the branch cut sweeps across the whole /c-plane, so that k moves from the lower to 
the upper sheet. Throughout the rest of the integration along 72, k remains on the 
upper sheet. 

Since W(z', k) is analytic in k away from the branch cut determined by z' , equation 



(2.7) defines <p(z, •) as an analytic function S z — > C whenever k does not lie on one 



of the dotted curves in Figure 2.1 associated with the movements of the branch cuts. 



Thus there are only two segments across which eft may have singularities: 

1. The branch cut [— iz, iz] of S z . The values of (f>(z, k~) to the left (right) of this 
cut glue together smoothly with the values of (ft(z, k + ) to the right (left). Thus, 
(ft is analytic as a map S z — > C near this cut. 

2. The contour T defined as the segment in the upper sheet of S z lying above 
[—ipo, ipo]. This segment, which arises from integration along the disk, is present 
only for the contour 72 (hence it lies only in the upper sheet). We will show in 



section 2.4 that (ft(z, •) does have a jump across F and we will express this jump 



in terms of the boundary values of U on the disk. 



Remark 2.1 1. The Lax pair (2.2 ) has singularities at the two branch points k = —iz 
and k = iz at which A = 00 and A = 0, respectively. However, as we verify in detail 
in appendix |A| the eigenfunction (j)(z, k) is still analytic near these points. 



2. The boundary condition (2.1b) ensures that the integration along the (-axis, 
which involves the degenerate Riemann surfaces, is compatible with the integration 
involving nondegenerate Riemann surfaces. 

2.4 The inverse problem 



Equations (2.7) provide the solution of the direct problem, namely they express (ft in 
terms of U for all k £ S z . In order to solve the inverse problem, we must find an 
alternative representation for (ft, namely we must express (ft in terms of an appropriate 
spectral function. This can be achieved by formulating a RH problem on S z (see 



Figure 2.2). It was shown in the previous subsection that 

(ft(z,k) is analytic for k £ S Z \T. (2.8a) 

Thus, in order to formulate a RH problem for (ft we must compute the 'jump' of (ft 
across F in terms of the boundary values of U on the disk. Introducing the notations 
(ft + and <j)~ for the values of (ft to the right and left of F, i.e., 

<ft + (z,k) :=<ft(z,(+0 + ik 2 ) + ) , <t>~(z,k) := (ft {z ,{-{) + ik 2 ) + ) , -p <k 2 <p , 
and denoting the jump by D(k), i.e., 

D(k) = (ft + (z, k) — (ft~(z, k), k€F, (2.8b) 



The case when k lies directly on the segment [—iz, iz] can be handled by slightly deforming the 
contour 71. 
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we will show that 

rPo -ik 2 (U+-U-)+p(U+-U^) 



D(k) 



k 2 



^2 



dp, k = ik 2 £ T, k 2 > 0, (2.9) 



where Up^ and denote the Dirichlet and Neumann boundary values of U on the 
disk, respectively, i.e. 



Uf :=U p (p±iO), 



U c (p±i0). 



In order to prove (2.9) we need to evaluate W(z,k) for k near the branch cut 
[C — ip,C + ip\- Suppose k = (C ± + ik 2 ) + , k 2 > 0, lies on the upper sheet. Then 
the expression (2.3) for A yields 

Mk=((+0+ik 2 ) + = 




Hence 



A 



+ A 



2k 2 



fc=(C+o+jfc 2 ) H 



k=(C+0+ik 2 )+ 



Similarly, 



+ A 



2ik 2 

2p 
2ip 



2k 2 



1 

A 



A 



k=((-0+ik 2 )+ 



k={(-0+ik 2 )+ 



-2ik 2 



-2ip 



P < k 2 , 
p > k 2 , 

P < k 2 , 
p > k 2 . 

P < k 2 , 
P > k 2 , 

P < k 2 , 
p > k 2 . 



(2.10a) 
(2.10b) 



(2.10c) 
(2.10d) 



Letting 

W R (z, k + ) := W (z, (C + + ik 2 ) + ) 
the definition (|2.5|) of W implies 



W L (z, k + ) := W (z, (C - + ik 2 ) + ) 



W R {z,k + ) 



1 

Vp^kl 



- [{ik 2 U p - pU ( ) dp + {ik 2 U c + P U P ) d(] , p < fc 2 , 
[(ifel/p - pi7 c ) dp + (i/c 2 £/f + pf^p) d(] , p> k 2 



(2.11) 



This notation turns out to be convenient and will be used frequently for functions of z and k: 
a superscript R or L on a function evaluated at a point (z, k + ) with k + lying on the branch cut 
— ip, f + ip] determined by z = p + iC,, means that the value of k + should be shifted infinitesimally 
to the right or left of the branch cut on the upper sheet before evaluation. 
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(2.12) 



and 

wH>X) = { wR{z R k+) : ' <fe ' 

1 ' \-W"(z,k+), p>k 2 . 

The jump D{k) of (j) across V can be determined by considering the integral J 2 W 
which defines (f)(z, k + ) for k = ±2e + ik2 and taking e — > 0. We choose the integration 
contour 72 so that when it passes along the disk, it lies a distance e above or below 
the disk (so that k remains on upper sheet throughout the integration along the disk) . 
The integrals along the subcontours of 72 which are not along the disk are continuous 
across T. Hence 



D(k) 



[0,Po] 



[W H (p - iO, k + ) - W H {p + iO, k+) 
W L (p - iO, k + ) + W L {p + iO, fc H 



(2.13) 



ker. 



In view of (2.12), this can be rewritten as 



D{k) =2 / [W R {p - iO, k^ 



W R (p + iO,k^ 



(2.14) 



Using in (2.14) the expression for W R given by (2.11), we find (2.9) 



The values of D(k) for k = ik?. € V, k% < 0, can be obtained from (2.9) by 
symmetry. Indeed, the relation 



X(z, k ) 



1 



implies that 



W(z, k^) 

It follows that <j) admits the symmetry 

<P(z,k ± ) = 

and so 

D{k) = D(k), 



\(z,k±) 
W(z,k±). 

k£T. 



(2.15) 
(2.16) 



Equations (2.8) constitute a RH problem for <j)(z, k) with jump across T. In order 



for the solution of this RH problem to be unique we need to supplement equations 



(2.8) with a normalization condition for cj)(z,k). The value of 4>{z,k + ) is defined by 



integration along 72 according to (2.7). The integration along the semicircle vanishes, 
so that 



-ioo, k + ) = 0. Therefore, in the limit k — * 00, we find 



lim cj)(z, k^ 

k— >oo 



lim 

k— >oo 



W(z,k^ 



dU = U{z). 



Similarly, 



lim (j)(z, k ) 

k^oo 



-U(z). 



(2.17) 
(2.18) 
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Figure 2.2 The contour T and the Riemann surface S z used for the RH problem. 



In particular, 



z, oo 



0,oo ). 



(2.19) 



Equation (2.19) provides the required normalization condition which ensures unique- 
ness of the solution of the RH problemj^] 



The unique solution of the RH problem (2.8) together with the normalization 



condition (2.19) is given by 



4>{z, k) 



1 

4ni 



D(k') 



X(z, k)(k + iz) 
X(z,k')(k' + iz) 



+ 1 



dk' 
k'-k' 



(2.20) 



Indeed, the right-hand side of (2.20) is an analytic function of (X,k) E S z for k ^ 



-ipo,ipo]. The standard Plemelj formulas imply that this function satisfies the jump 



condition (2.8b) across T in the upper sheet, whereas it does not jump on the lower 



sheet. Finally, the condition (2.19) is a consequence of the relations X(z,oo + ) = 
— X(z, oo~) = 1. 

Equations (2.7) express 4> in terms of U (the solution of the direct problem), 
whereas equation (2.20) expresses <ft in terms of D{k) (the solution of the inverse 
problem). Using these two different representations of <p it is straightforward to 



compute U in terms of D{k): Substituting the representation (2.20) into (2.17), we 
find 

U(z) = ±-l -, ^ . r dk. (2.21) 
v ' 4m J T X(z,k)(k + iz) v ' 



This can be written as 



U(z) 



1 

4iri 



D{k) 



-dk, 



where the branch with positive real part is chosen for the square root. 



(2.22) 



Remark 2.2 In the above discussion we assumed £ > 0. The case of £ < is similar. 



If 4>\{z,k) and <$>i{z, k) are two solutions of (2 



satisfying (2.191, then cj>i 
everywhere on 5 Z . Hence <f>\ — (j>2 is a constant. Since (<f>i — 4>2){z, oo + ) = —{4>i — 
constant is zero. 



— <f>2 is analytic 
4>2){z, oo~), this 
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2.5 The spectral functions 



The expression (2.9) for D(k) involves both the Dirichlet and the Neumann boundary 
values. However, for a well-posed problem only one of these boundary values is 



specified. In sections |2.6| and 2.7 we will use equatorial symmetry and the global 
relation to express D{k) in terms of either the Dirichlet or the Neumann boundary 
values. In this connection, we introduce the following definitions: For k = ik 2 with 
< k 2 < POj the spectral functions R{k) and S(k) are defined by 



and 



R(k) :-- 



S(k) :-- 



W R {p + iO,k + ) 



[0,k 2 



W R (p + iO,k + ) 



k2 ik 2 U+ - pU. 



po ik 2 U+ 



-p- 



dp 



dp. 



(2.23a) 



(2.23b) 



[k2,po] J hi 

Note that the functions Re(i?) and Re(S') are defined in terms of the Dirichlet and 
Neumann boundary values of U on the disk, respectively. 

2.6 Equatorial symmetry 

For a variety of physically significant BVPs, we expect the spacetime metric to be 
symmetric with respect to the equatorial plane. In terms of the Ernst potential this 
means that f(z) = f(z), so that in the static case U(z) = U(z). 

Proposition 2.3 Assume that U{z) is equatorially symmetric, i.e. U{z) = U{z). 



Then the spectral functions D(k) and S(k), defined by (2.9) and (2.23b), are related 
by 



D(k) = -4Re(S(k)). 
Furthermore, for z on the disk, z = p ± iO, 



W L (p - iO, k + ) = W R {p + iO, k+), 



keT. 



Proof. The symmetry U(z) = U{z) implies that 

U p (z) = U p (z), U c (z) = -U c (z). 

In particular, the boundary values above and below the disk are related by U~ 
and U, 



(2.24) 

(2.25) 

(2.26) 
= U~ 



with (2.23b), we find (2.24). 



Uf . Using these relations in (2.9) and comparing the resulting equation 



Let W = Widp + W 2 d(. By the definition (|2.5|) of W, we have 



I ■ A j U„{z)+ir x - \ )l 



\=X(z,k+) 



Using (2.26) together with the identity 

1 

~ \{z,k+ 



X(z,-k ), 



(2.27) 
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we find 



A [) U p (z)-i(-\+ l x )l 



J X=X(z,-k-) 

But the right-hand side of this equation is precisely —W\(z, —k~). The relation 

\( Zj k + ) = -X(z,k~), 

implies that 



W{z,k + ) = -W(z,k~). 
Thus, utilizing the symmetries ( |2.15[ ) and (2.28), we find 



(2.28) 



Wi(z, k + ) = -Wi(z, -k~) = Wi(z, -k + ) = Wi(z, -k+). 



The identity (2.25 ) now follows by taking z = p+iO and observing that if k + 



lies just to the left of T, then —k + = + ik2 lies just to the right of T. 



-0+ik 2 
□ 



Equations (2.23b) and (2.24) express the jump D(k) in terms of the Neumann 



boundary values. Thus, inserting equation (2.24) into the right-hand side of (2.22) 



we immediately obtain the solution of an equatorially symmetric Neumann boundary 
value problem. In order to express D{k) in terms of the Dirichlet boundary values, 
we will use the global relation. 

2.7 The global relation 

The global relation is an algebraic equation satisfied by the spectral functions. It 
expresses the fact that the boundary values are related and cannot be independently 
prescribed. For our problem the global relation can be derived from the equation 



W = 0, 



(2.29) 
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where 7 is the contour encircling the physical space depicted in Figure 2.3 The 



identity (2.29) is obtained by applying Stokes' theorem together with the fact that 



dW = in view of (2.4). 



Since W depends on k, (2.29) is a family of relations parametrized by k. We 



choose k = e + iki and specify the integration contour so that when it passes along 
the disk, it lies a distance 2e above or below the disk (so that k changes sheets during 



the integration along the disk). Equation (2.29) then reads 



/ W(z,k + )+ [ W R (p-iO,k + ) (2.30) 

J[-ioo,0] J[°>Po] 

+ [ W L (p + iO,k-)+ [ W(z,k~) = 0. 

J[ P 0,0] J[0,ioo] 



13 



Im k 




Re k 



Figure 2.3 The integration contour encircling all of physical space used for the global 
relation and the associated movement of the endpoints of the branch cut. 



We claim that the two integrals along the C-axis can be simplified as follows: 



[-ioo.O] 



W(z,k + ) = U(-iO), 



[0,ioo] 



W(z,k~) = U(+iO). 



(2.31) 
-1) for 



Indeed, for z = i( the form of (2.4) is particularly simple since A = 1 (A 
all k on the upper (lower) sheet. We find, for any £ £ R and k G C, 

W(i(, k + ) = dU, W(i(, k~) = -dU. 

Together with the initial conditions 4>(ioo, k~) = (p(—ioo, k + ) = 0, this shows (2.31 ). 
This yields the following expression for (2.29): 

U(-iO)+ I W R (p-iO,k + )+ [ W L (p + iO,k-) + U(+iO) = 0. (2.32) 

In the presence of the equatorial symmetry, equation (2.32) can be simplified 
further. Indeed, in view of Proposition 2.3 and the symmetries (2.12) and (2.28), we 
have 



W R (p - iO, k + ) 



[0,po] 



[0,fc 2 ] 



W L (p-iO,k + ) 



[ka,Po] 



I 

Apo,o] 



W L {p - iO, k + ) = R{k) - S(k), 
(2.33) 

W L (p + iO,k-)= I W R {p + iO,k + )- [ W R {p + iO,k + ) = -S{k) + R(k). 

^0,0] J\paM J[kn,Q] 

(2.34) 

Substituting these expressions into (2.32), we find the following result. 

Proposition 2.4 Suppose that U{z) obeys the equatorial symmetry U(z) = U(z). 
Then the spectral functions R(k) and S(k) defined by (2.25ty satisfy the global relation 



U{+iO) + Re(R(k)) = Re{S{k)). 



(2.35) 



14 



Since Re(i?) and Re(S') involve only the Dirichlet and Neumann boundary values 
of U on the disk, equation (2.35 ) defines the Dirichlet to Neumann map. In particular, 

15) yield 

D{k) = —4 (U (+i0) + Re(R(k))) , (2.36) 



equations (2.24) and (2.35) yield 



which expresses D in terms of only the Dirichlet boundary values. 

We summarize our discussion of the linear problem in the following theorem. 



Theorem 2.5 Let U{p,C) be a real-valued solution of the axisymmetric Laplace equa- 



tion \2. 1 o| ) in the exterior disk domain T> satisfying the boundary conditions (2.1b) 
and ( fg.icj ). Suppose that U is equatorially symmetric, U(p,() 
U(p, Q admits the integral representation 



U(p,-C). Then 



U(p,Q 



1 

4iri 



■i P» 



D(k) 



V(k-C) 2 + P : 



-.dk. 



(2.37) 



where the branch with positive real part is chosen for the square root, and the function 
D(k) is given in terms of the Dirichlet and Neumann boundary values of U on the 
disk by the following expressions respectively: 



D(k) 
and 



-4 U{+iO) + \k 2 \ 



^ U p (p + iO) 



dp 



k = ik 2 , \k 2 \ < po (2.38) 



D(k) 



p° pU c {p + iO) 
Ifcal V? 



dp, 



k = iko 



\k 2 \ < Po- 



(2.39) 



It can be verified directly, using Abel transforms, that the integral representation 



(2.37) indeed yields the correct boundary values for U(p,(), see appendix |B| 



3 The Ernst equation 
3.1 Lax pair 



The elliptic Ernst equation (1.1) admits the Lax pair 

U z {z,k) = U{z,k)<S>{z,k), 
\<5> s {z,k) = V{z,k)<S>{z,k), 



(3.1) 



where z = p + iC, the function $(z, k) is a 2 x 2-matrix valued eigenfunction, and the 
2 x 2-matrix valued functions U and V are defined as follows: 



U 



B XB 
XA A 



A 



fz 

/ + /' 



V 



B 



A 



B 



fz 

f+r 
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with A = \(z, k) given by (2.3 ). For each z, &(z, •) is a map from the Riemann surface 
S z to the space of 2 x 2 matrices. As before, we use the notation $(2, k) to denote 
$(p,C,k), etc. 



We can write the Lax pair (3.1 ) as 

d$ = W$ 

where W is the one-form 
W =Udz + Vdz 

1 ( u 



(3.2) 
(3.3) 



/ + / \h [(i + A)/p + <Q-A)/ c ] 
1 



1 ry 1 
2 



+ 



[a + A)/p + <(i-A)/ c ]\ 

/p / 

H(x + A)/c-*( A -A)/ P ] 



/ + / Vf [(x + A)/ C -i(i-A)/„] /c 

For a 2 x 2-matrix A, we let [A]i and [A] 2 denote the first and second columns of 
A, respectively. We define an eigenfunction 3>(z, fe ) as the solution of (3.2) which 
satisfies the initial conditions 



lim [$(z, k ) 
lim [<f>(z,k + ) 



-1 



for all k on the lower sheet, (3.4a) 
for all k + on the upper sheet. (3.4b) 



These initial conditions are convenient because they lead to the following symmetry 
properties of 3>: 

$(z,k + ) = <r 3 $(;z,fc-)<7i, $(z,k + ) = <ji$(z, k+)a 3 . (3.5) 



Indeed, these symmetries are a consequence of (3.4), as well as of the following non- 
linear analogs of equations (2.15) and ( |2.28 ): 



W{z, k + ) = a 3 W(z, k~)a 3 , W(z, k + ) = a x W{z, k+)a 1 . 



(3.6) 



Remark 3.1 The Lax pair (3.1) may have singularities at points where Re/ 



0. Physically these points make up the ergospheres of the spacetime (within these 
surfaces there can be no static observer with respect to infinity). We will henceforth 
assume that no ergospheres are present, although we expect our discussion to apply 
to many cases of physical interest as long as Re / > on the (-axis. 

3.2 The direct problem 

Assume that z = p + iQ with £ > 0. Just like for the linear problem, we can express 
5>(z, k) for all k on S z by using integration with respect to the contours 71 and 72 in 
Figure [2~7[| i.e. 



*(*,AT)= lim $(z,k~)+ [ (W®)(z',k), 
$(z,k + )= lim $(z,k~)+ [ (W$)(z',k). 



(3.7) 
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For the linear problem there was a jump across T = [—ipo,ipo] in the upper sheet. 
For the nonlinear problem we will find that jumps both across T + and T~ , where 
r + and r~ denote the coverings of T = [—ipo,ipo] in the upper and lower sheets, 
respectively. 



The value of [&(z, fc~)]i as z —* zoo is fixed by (3.4a), thus an argument similar to 
that used for the linear problem applied to the contour 71 shows that [$] 1 is analytic 
on the lower sheet. A similar argument applied to the contour 72 shows that [$]i is 
analytic on the upper sheet away from the contour r + . From the first symmetry in 



(3.5), it follows that [$]2 is analytic on the upper sheet and on the lower sheet away 



from r . 

3.3 The inverse problem 

In order to formulate a RH problem for we need to determine the 'jump matrices' 
across T + and T~ . Suppose that we can express the jump matrix D, defined by 

$-(z,k) = $ + (z,k)D(k), k£T + , 



in terms of the boundary values of / on the disk. Then, the first symmetry in (3.5) 
immediately gives the following expression for the jump across 

9~(z,k) = <^ + (z,k)a 1 D(k)a 1 , k G V. 

In the linear problem D was computed directly in terms of the relevant boundary 
values. However, for the nonlinear problem it is more convenient to proceed in three 
steps. In step 1 we express the values of $ on the £-axis in terms of two spectral 
functions F(k) and G(k). In step 2 we use the contour 71 to express the jump D{k) 
in terms of F(k) and G(k). In step 3 we use the contour 72 to relate F and G to the 
boundary values of / on the disk. Combining steps 2 and 3 we obtain an expression 
for D{k) in terms of the boundary values of /. 

Step 1. The values of on the (,-axis can be expressed in terms of two spectral 
functions F{k) and G{k) as 

=(«g C>0. tEC , (3.8a) 

•(*,*-)= (J _^,)(J <>0. k^C. (3.8b) 

•«,*-)= (j %){gw !)■ c<0 - tec (3 - 8d) 

The functions F(k) and G{k) have the following properties: 

• F and G are unique functions of k G C, i.e. viewed as functions on S z they 
satisfy 

F{k + ) = F(k-), G{k + ) = G(k~), k e C. (3.9) 
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F(k) and G(k) are analytic for k 6 C \ 
Under the conjugation k >—^k, F and G obey the symmetries 
F(k) = G(k) = -G(k), keC. 



(3.10) 



In the limit k — ► oo, 

F(k) = 1 + 0{l/k), G{k) = 0{l/k), 



k — > oo. 



(3.11) 



Proof. For z = i£, A = 1 for all k + on the upper sheet and A 
the lower sheet. Thus 

1 fdf dp 
f + f\df df / 



W(i(,k + ) 



-1 for all k on 



(3.12) 



We infer that two independent vector-valued solutions of d<£ = W(i^, k + )& on the 
C-axis are 



/(*) 
f(iC) 



(3.13) 



These solutions together with the initial condition (3.4b) imply that equation (3.8a) 



holds for some functions F(k) and G(k). Equation (3.8b) follows from the first 



symmetry in (3.5). The value of $ at z = — ioo is obtained from the value z = ioo by 



integrating along a large semicircle at infinity. Thus, using (3.8b) and the fact 
that vanishes for large z, we find 



lim $(z, k A 



lim $>(z,k 



i n i G(k) 

1 -lMo F(k) 



(3.14) 



Using this initial condition together with (3.13), we find (3.8c). Finally, (3.8d) follows 
by symmetry. 



By definition, F and G satisfy (3.9). The analyticity properties of F and G follow 



from the analyticity properties of (f)(z,k). The symmetry properties (3.10) are a 



consequence of (3.8) and of the second symmetry in (3.5). In order to prove (3.11) 



we note that in the limit k + — > oo, equation (3.12) holds for all z, i.e 



Thus, 



lim W(z, k A 



lim $(z, k + ) 

k^oo 



1 fdf df 



f + f\df df 



-i 



a 
b 1 



4 This property holds because of our assumption that Re / > on the £-axis (the solitonless case). 
More generally, F and G will have poles which generate solitons. For example, for the Kerr black hole 
(which is a two-soliton solution), F(k) and G(k) have two poles each on the real axis corresponding 
to the fact that the real part of / vanishes at two different points on the ("-axis. 
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for some constants a, b 6 C. Evaluating this equation at z = —ioo and comparing 
with (3.14), we find that a = 1 and 6 = 0. Thus, 

$(z,k + ) 



/(*: 
/(*: 



+ 0(l/fc), A;^oo. 



Comparing this equation with (3.8a), we find (3.11). 



□ 



The functions F{k) and G{k) jump across T = [—ipo,ipo]- Let F + ,G + and 
F~,G~ denote the values of F and G for k to the right and left of T, respectively. 
These jumps are related to D, as shown below. 



Step 2. The jump D across T + is given in terms of F and by 



F-(lfe) 0\ 
G-(k) I]' 



(3.15) 



Proof. Integration along the contour 71 with k on the upper sheet gives 

(W$)(z',k+), 



•<**♦>-(! -OP " l + 



71 



where we have used (3.8a) to determine the initial condition. Since the branch cut 
never comes near T when integrating along 71, the function W in the integrand is 
continuous across T. Thus the jump A<I> := <I> + — satisfies 



A$(z,fc+) 



1 1 \ / (F+(k) 
1 -I [\G+(k) 1 



F-(k) 
G-(k) 1 



+ 



(T^A$)(z',A: + ). 



71 



This shows that A<I> satisfies the same differential equation as $ + , so that there exists 
a 2 x 2-matrix valued function D{k) independent of z such that 



A$(z,k + ) = $ + (z,k + )(I - D(k)) i.e. $~ = *+D(fc). 



(3.16) 



Evaluation of (3.16) at z — > ioo, yields the following equation: 



1 1 \ ( (F + {k) 
1 -lJ^G+tA;) 1 



F-(Jfe) 
G-(Jfe) 1 



1 1 \ /F+(A;) 



1 -1 



U + ; fc ) ij( j -^))- 



Solving this equation for D, we find (3.15). 



□ 



Remark 3.2 Equation (3.15) shows that the second column of D{k) is given by 
(0, 1) . This is consistent with the observation made earlier that the second column 
of $ has no jump across T + . 
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Step 3. The functions and satisfy, for k = i/c2, < k2 < po, the following 
equation: 



1 G~(k) 
F-(k) 



1 G + (k)\ (F+{k) 0\ 1 (F~[k) 0\ _ 1 (I 1 



C(k), 



v F+{k)J \G + {k) lj \G~(k) I J 2 VI 
where the 2 x 2-matrix valued function C{k) is defined b}J^] 
C(fc) = f [({W R - W L )<S>- L ) (p - iO, k + ) - ({W R - W L )<S>- L ) (p + iO, k + )] 



(3.17) 



(3.18) 



Proof. Consider the expression for $(z, k + ) given in (3.7) by integration along 72. 
Let k + = ±2e + iki and choose the integration contour 72 so that when it passes 
along the disk, it lies a distance e above or below the disk (so that k remains on the 
upper sheet throughout the integration along the disk). Let j' 2 denote the part of 
72 which does not run along the disk. Then the W in the integrand is continuous 
across T + during the integration along 7 2 . It follows that A<£ satisfies the following 
equation: 

A$(z,k + ) = lim A$(z,k~) + f (WA&)(z', k) 

z->ioo Jy 2 

+ f [(W$) +R (p-iO,k + ) - (W$) +R (p + iO,k + ) 
J[0,Po] 

- (W$y L {p - iO, k + ) + {W$>)~ L {p + iO, k + )] . 
The integral along the disk can be written as 

C{k)+ [ [(W R A$)(p-iO,k + ) - (W R A$)(p + iO,k + )], 
J[0,p ] 



where C(k) is given by the right-hand side of (3.18) with the lower integration limit 
replaced with 0. Since W L = W R for p < &2, we see that C(k) can be expressed as 
in (3.18). We infer that A3> and <£ + satisfy the similar equations 

A$(z,k + )= lim A*(*, *T) + C(fc) + [ (WA$)(z', k), keT 



and 



$+(z,k + )= lim $ + (z,k + )+ f (W$ + ){z',k), keT, 



where W is W R for the integrations along the disk. Comparing these equations and 
using (3.8b) and ( |3.16 ), we find 



1 1 



1 G+{k)\ (I G~(k) 



+C(k) 



1 1\/1 G+(k) 
1 -lMo F + (k) 



(I-D(k + )). 



1 -1) v v° F+ i k )J v° F ( k ) 

The superscript — L indicates that the function should be evaluated with k lying on the — side 
of r and to the left of the branch cut. 
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Using the expression (3.15) for D and simplifying, we find (3.17) 



□ 



Remark 3.3 1. The matrix C(k) in (3.18) can be computed in terms of the boundary 



values of / on the disk via only linear operations. 



2. In (3.18) we have assumed that k 2 > 0; the jump for — po < k2 < can be 



obtained via the symmetry D{k) = a^D{k)az- 



3. Equations (3.17) and (3.18) show that, in the case that W is free of singularities, 



there is no jump at the endpoints of T, i.e. D(±ipo) = I. 



In order to compute C(k) we need the nonlinear analogs of equations (2.11 ) and 



(2.12 ), which express W in the vicinity of the branch cut. On the disk, these equations 



arc 



W R {p±iO,k- 



disk 



fp 

f+f I ihfp-pft 



2_ n 2 
2 P 



_ fp 

f+f \ ik2j P -Ph 



ik2f P ~pf( " 
iy/kl-p 2 

fp 

fp 



dp, 



z=p±i0 



dp, 



p < k 2 



P> k 2 



(3.19) 



z=p±i0 



and 



W L (p±iO,k^ 



(w R ( P ±iO,k+), P <k 2 , 
\a 3 W R (p±iO,k+)a 3 , p>k 2 . 

The above results are summarized in the following proposition. 



(3.20) 



Proposition 3.4 Let f(p,C) satisfy the Ernst equation (1.1) in the exterior disk 



domain T>. Suppose that f is asymptotically flat, i.e. f{p,C) ~ ¥ 1 as P ' + C —> oo, 
regular on the rotation axis, i.e. /p(+0, £) = for all £ ^ 0, and that no ergospheres 
are present, i.e. Re/ > in T>. Then the solution f(p,C) can be expressed in terms 
of both the Dirichlet and the Neumann boundary values on the disk as follows: 

1. Use the boundary values of f to find W R {p ± iO, k + ) and W L (p ± iO, k + ) from 



equations (3.19) and (3.20). 



2. Find $ L (p ± iO, k + ) for < p < po and < k 2 < po by solving the ordinary 
differential equation 



<5>- p L {p ± iO, k + ) = (W<5>)~ L (p ± *0, A; H 



< p < Po, 



together with the initial conditions 
[<I>- L (-iO,A: + )] 1 : 



/Ho) 
/Ho) 



[$- L (+iO,k + )] 2 



as well as the following continuity condition at the tip of the disk: 
$- L (p - iO, k + ) = §- L { Po + iO, k + ). 
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3. Use the expressions for W R , W L , and $ L to compute C(k) according to (3.18). 



4- Find D(k) from (3.15) and (3.11). The entries of D are rational functions of 
the entries of C . 

5. Compute $ in terms of the spectral function D{k) via the solution of the fol- 
lowing RH problem: 

• $>(z, k) is an analytic function of k £ S z away from k £ r+ . 

• Across r + , $ satisfies the jump condition 

*-(z, k) = $ + (z, k)D(k), k £ F + . 

• Across r~ , $ satisfies the jump condition 

$~(z,fe) = ^ + {z,k)aiD(k)a 1 , k€T~. 

• As k — > 00, $ satisfies 



lim [$(z, fc" 

fc— >oo 



lim [$(z, k + )} 2 



• $ obeys the symmetries 

$(z,k + ) = <r 3 ^(z,k~)<n, 
6. Find f{z) from the equation 



<fr(z, k + ) = cti®(z, k+)a 3 . 



f{z) = lim (*(*,fc+)) 2 i. 



(3.21) 



(3.22) 



(3.23) 



Remark 3.5 The solution of the RH problem in Proposition 3.4 is unique, despite 
the nonstandard form of the normalization condition (3.21 ). Indeed, if $1 and $2 are 
two solutions of this RH problem, then the functions $j, j = 1,2, defined by 

$j(z, k) = §j(z, oo + )- 1 ^ j {z, k), j = 1, 2, 

satisfy the same jump condition as $1 and $2 together with the standard normaliza- 
tion condition $j(z, oo + ) = /. Thus $1 = $2- The equality $1(2, oo~) = $2(2, oo~) 
together with the first symmetry in (3.22) yield 



$l(z,oo+) 1 0-3^1(2;, oo + )cji = $2(^,00+) V 3 $2(2,oo + )cji 



Equation (|3.21|) together with second symmetry in (|3.22|) imply that 

1 



-1 



i = 1,2, 



(3.24) 



(3.25) 



where f\ and f 2 are complex- valued functions. Equations (3.24) and (3.25) imply 
that $1(2,00+) = $2(2,00+). Hence 



$1(2, k) = $i(z, oo + )$i(z, k) = $2(2, oo + )$ 2 (2, k) = $2(2, k). 
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4 Equatorial symmetry and the global relation 



Proposition |3.4| expresses the solution of the Ernst equation in the domain T> via the 
solution of a RH problem formulated in terms of both the Dirichlet and Neumann 
boundary values of / on the disk. Since only a subset of the boundary values can be 
specified for a well-posed problem, the solution formula is not yet effective. Recall 
that for the linear problem, by using equatorial symmetry together with the global 
relation, we were able to determine the jump data for the RH problem in terms of 
only either the Dirichlet or the Neumann boundary values. Similarly, for the non- 
linear problem, for equatorially symmetric solutions whose boundary values possess 
a sufficient amount of symmetry (such boundary values are called linearizable) , the 
unknown boundary values can be eliminated. 

In this section we analyze the consequences of / being equatorially symmetric, 
i.e. f(z) = f(z), and derive for this case the global relation satisfied by the spectral 
functions. The linearizable boundary conditions are analyzed in section [5] 

4.1 Equatorial symmetry 



The nonlinear analog of Proposition 2.3 is the following. 



Proposition 4.1 Assume that f is equatorially symmetric, i.e. f(z) = f(z). Then, 
for k G T, 

W L (p - iO, k + ) = W R (p + iO, k+j. (4.1) 

In particular, there exists a 2 x 2-matrix valued function K{k), k £ V, independent 
of p such that 



$- L (p-iO,k + ) = $>+ R (p + iO,k+)K(k), keT. (4.2) 



The spectral functions F(k) and G(k), introduced in section 3.3, are related to K{k) 
as follows: 

K{k) = A^fk)a 1 A-{k)a 1 , k£T, (4.3) 

where 



K G(k) 1 

and A± denote the values of A to the right and left of T respectively 



Proof. Let W = W\dp + W2d(. Then, for any z, the definition (3.3) of W yields 



o-i<t 3 Wi(z, k + )a 3 a 1 = — — l -==a 1 a 3 
f{z) + f{z) 

\\ [(i+A)/„(z) + i(£-A)/ c (f)] " f p {z) 



\ <73<7i. 

A=A(z,fc+) 
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Using the equatorial symmetry, the right-hand side of this equation can be written 

as 



1 



/(*) + f(z) VJ [-(i + X)f p (z) + i(i - A)/ C (z)] 



-(i + A)/ p (z) + i(i-A)/ c (z) 



In view of (2.27), this expression equals W\(z,—k ). Therefore, utilizing the two 
symmetries in (3.6), we find 



ax<r 3 Wi(z, k + )a 3 a! = Wi(z, -k ) = a 3 Wi(z, -k + )a 3 = a 3 axWi(z, -k+)aia 3 . 

This proves that 

Wi{z,k + ) = Wx(z,-k+). 



Evaluation of this identity for z = p + iO on the disk yields (4.1). Indeed, if k + = 
— + ik2 lies just to the left of T, then —k + = + iki lies just to the right of T. 

It follows from (4.1) that <&~ L (p — iO,k + ) and $> +R (p + iO, k + ) satisfy the same 
differential equation. This establishes the existence of K{k). 



Evaluating (4.2) at p = 0, we find 



K{k) = <S>+{+iO,k+) 1 $-(-i0,k + ). 



(4.5) 



In view of the axis values (3.8) of and the equatorial symmetry, equation ( |4.5| ) 
yields (Ok. □ 



A=A(z,fc+) 



4.2 The global relation 

In this subsection we derive the global relation satisfied by the spectral functions 
F(k) and G(k). 



Proposition 4.2 Assume that f(z) = f(z) and define K(k), k 6 V, by (4-2). Then 
K(k) satisfies 

axKik) = K{k)a u k&T. (4.6) 

Proof. For any z, it holds that <&(z, k + ) = &(z, k~) for k = iz since the two sheets 
of the Riemann surface coincide at this branch point. Consequently, in view of the 
first symmetry in (3.5), 



$(z,k + ) = az<S>(z,k + )o-i, 



iz. 



(4.7) 



Evaluating (4.7) for z on the disk, we are able to establish (4.6). Indeed, choosing z = 
ki + i0 and z = k% — iO in (4.7), and introducing the notation <&k 2 := ^ + (k 2 -MO, ik^ ), 
we find 

$fc 2 = CT 3 $ fc2 <Ti (4.8) 



and 



$ (k 2 - i0,ik2 ) = a-z® (&2 - iQ,ik£)o-i, 



(4.9) 
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respectively. Moreover, evaluating (4.2) at p = kz, we find 

* - (A& -iO,k + ) = $> k2 K(k), ker. 



Equations (4.9) and (4.10) yield 



$ k2 K = az$k 2 Kai. 



(4.10) 



(4.11) 



Using (4.8) to repla ce o~3<&k 2 with $fe 2 <Ji, and then premultiplying both sides by 
CTi^r 1 , we find (O. 



□ 



Propositions 4.1 and |4.2| imply the following result. 



Proposition 4.3 Suppose that f is equatorially symmetric, i.e. f(z) = f(z). Then 
the spectral functions F(k) and G(k) satisfy the following relation, which will be 
referred to as the global relation: 



A + {k)a 1 A~ 1 {k)a 1 = a 1 A^(k)a 1 AZ 1 (k), k£T, 



(4.12) 



where A(k) is defined in terms of F(k) and G(k) by equation (4-4) '• 



Remark 4.4 The global relation (4.12) can also be derived by using the nonlinear 
analog of equation (2.29), but the derivation presented here is simpler. 



5 Linearizable boundary conditions 



In Proposition 4.3 we derived, under the assumption of equatorial symmetry, the 



global relation (4.12) satisfied by the spectral functions F(k) and G(k). In this 
section, we will show that if the corotating potential /' = /q is constant on the disk, 
then the spectral functions F{k) and G{k) satisfy an additional important algebraic 
relation. Furthermore, we will show that these two algebraic relations satisfied by 
F{k) and G{k) yield an auxiliary RH problem for F and G with jump data given in 
terms of only the known boundary value /q. Moreover, this matrix RH problem is 
diagonalizable and reduces to a scalar RH problem, which can be solved explicitly. 
In this way we recover the celebrated Neugebauer-Meinel disk solutions [T8]-|20 . 

Before considering the consequences of the boundary condition /' = /q, /q con- 
stant, on the disk, we first recall the physical origin of the Ernst equation and describe 
the corotating potential /' (see [21] for further details). 

5.1 The Ernst equation 

In canonical Weyl coordinates the exterior gravitational field of a stationarily rotating 
axisymmetric body is described by the line element 

ds 2 = e- 2U (e 2K (dp 2 + d( 2 ) + p 2 dp 2 ) - e 2U (dt + ad^) 2 , (5.1) 
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where p, £, ip are cylindrical coordinates, t is the coordinate time, and the metric 
functions U, K, a depend only on p and (. In these coordinates, the Einstein field 
equations consist of the equations 



u pp + u (C + -u p 



~2f K + a c)> 



(P 



0, 



(5.2) 
(5.3) 



together with two equations involving k. The condition that the metric is regular at 
the rotation axis implies that 



0, k — ► as p — ► 0, 



(5.4) 



whereas the condititon that the line element approaches the Minkowski metric at 
infinity (asymptotic flatness) implies that 



U -» 0, a 



0, k — ► as /o 2 + C 2 



OO. 



In view of (5.3), it is possible to introduce a function b(p, () such that 

-4Ui 



-pe~ 4U b p , 



and 



(p e - 4 %) p + (pe- 4 ^ c ) c = 



4(7, 



(5.5) 



(5.6) 



(5.7) 



Letting / = e 2U + ib, equations (|5.2|) and (|5.7|) combine into the single Ernst equation 



( 1.1 ). Moreover, as a consequence of the Ernst equation, the compatibility conditions 

Ken ar e automatically satisfied and the metric functions a 



a^ p and n 



and k can be determined by integration of the two equations (5.6) and the two 



field equations for k, respectively. Thus, the vacuum Einstein field equations in the 
stationary axisymmetric case are equivalent to the Ernst equation. 

5.2 Corotating coordinates 

Let us introduce the corotating coordinates (p' , (' , <p' ,t') bjj^] 

9=9, (' = C, ¥>' = <P ~ tit, t' = t, 
where Q is the constant angular velocity of the body. In these new coordinates, the 



metric (5.1 ) retains its form and the corotating metric functions U', a' , k' are related 

(5.8) 



to U, a, k via 



e 2U ' = e 



(l-na')e 



2„2„-4U 



[(1 + na) 2 - SFffe 



k' -U' = k-U. 



Since the form of the line element is invariant, the field equations retain their form 
in the corotating system. Thus, we may introduce a corotating Ernst potential /' by 



3 We will use primes to denote corotating quantities. 
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f = e + ib' and the Ernst equation retains its form in the corotating system as 
well. 

The Lax pair in the corotating system involves an eigenfunction and the one- 



form W' defined by replacing / with /' in (3.3). The eigenfunction <J>' is defined as 



the solution of d& = W'$>' which satisfies the initial conditions (3.4) with replaced 
by <£'. It can be verified 



where 



that the corotating eigenfunction $' is related to $ by 
&(z,k) = A(z,k)$(z,k), k£S z , (5.9) 



A(z, k) = (1 + na)I - npe~ 2V a 3 + i(k + iz)Vte~' 2U (-<t 3 + X(z, k)(na 3 ). (5.10) 



2Ui 



Thanks to (5.5 
conditions ( 



the relation (5.9) is consistent with the requirement that the initial 



3.4) should retain their form for <J>'. 



The spectral analysis of the corotating Lax pair is similar to that of (3.1 ), except 



that /' and / satisfy different boundary conditions. For example, equation (5.8) 
implies that Re/' ~ 1 — J7 2 p 2 as p — » 00, reflecting the fact that the metric is no 
longer asymptotically flat in the corotating system. Thus, a given BVP may possess 
additional symmetries in one of the two coordinate systems. In the next subsection 
we investigate the consequences of /' being constant along the disk (this corresponds 
to more complicated boundary values for /). 

5.3 /' constant on the disk 

We consider the condition /' = /q, f' Q constant, on the disk. It turns out that 
the condition of /' = on the disk implies an important relation (see Proposition 



5.3) satisfied by the spectral functions F(k) and G(k). The following analysis is 



conceptually similar to the analysis presented in section [4| Propositions 5.1 5 .3| below 
are the direct analogs of Propositions 4.1j|4.3 The resulting algebraic relation is 
however independent of the global relation (4.12). In this subsection we will not 
assume that / is equatorially symmetric. 

Proposition 5.1 Assume that f = on the disk. Then, for k £f, 



W' R (p + iO, k + ) = a z axW ,L {p + »0, /c+)ctict 3 . (5.11) 

In particular, there exists a 2 x 2-matrix valued function Q{k), k € V, independent of 
p such that 

<S>'+ L {p + iO,k+) = a 3 a 1 ^' +R (p + iO,k + )Q(k), keT. (5.12) 



(5.13) 



The spectral functions F(k) and G(k) introduced in section 3.3 are related to Q(k) as 
follows: 

Q{k) = A~ 1 {k)B- 1 A- 1 {k)a x a 3 A{k)BA + (k), keT, 
where A+ is defined in terms of F and G in Proposition 4-1 an^ 

7(+io) i\ 0.//H0) 1 

J(+i0) -1)' • \f(-i0) -1. 



D 



A(k) := A(+i0, k + ). (5.14) 



7 In the presence of equatorial symmetry B will be the complex conjugate of B. 
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Proof. By the second symmetry in (3.6), 



o x W ,L (p + iO, k+)a 1 = W' L (p + iO, k + ). 



Thus, for z = p + iO on the disk, the corotating analog of the definition (3.3) of W 
shows that the right-hand side of (5.11) is given by 

a 3 W /L (p + iO,k + )a 3 (5.15) 



1 



/' + /' 



f 

(x + A)/;+*(i-A)/; 



(i + A)/;+^-A)/^ 
f 



dp. 



On the other hand, 
W' R {p + iO,h + ) 
1 



f'o 



f> + f> ^1 ( 1 + A)/ / + i( 1 _ A)/ / 
For any z, 

A(z,fc+) 

For z = p + iO, this equation yields 

X R (p + iO, k + ) 



f 



z=p+iO 
X=X L (z,k+) 



(5.16) 
dp. 



z=p+iO 
X=X R (z,k+) 



\(z,(-k + 2()- 



X L (p + iO,k+)' 

Using this relation to replace X L (z,k + ) with l/X R (z,k + ) on the right-hand side of 
(5.15), and then subtracting the resulting equation from equation (5.16), we find 

W' R {p + iO, k + ) - a 3 a 1 W' L (p + iO,k+)a 1 a 3 

dp. 



1 







/'+/' V(x + a)/; 



:i + A)/; 







z=p+iO 
X=X R (z,k+) 



Setting /' = in this equation, we find (5.11 ). 



It follows from (5.11 ) that the functions a 3 ai& +R (p+iO, k + ) and $ / + L (p + iO, k+) 
satisfy the same differential equation. This establishes the existence of Q{k). 



Evaluating (5.12 ) at p = 0, we find 



Q{k) = & + (+i0,k + )- 1 o- 1 o- 3 $'+(+i0,k+). 



(5.17) 



In view of the axis values (3.8) of $ and the definition (5.9) of A, equation ( 5.17| ) 
yields (15131). □ 
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Proposition 5.2 Assume that f' p = 0on the disk and define Q(k), k ET, by (5.12). 
Then Q(k) satisfies 

axQ{k) = -Q(k)ai, k^T. (5.18) 



Proof. Evaluating (5.12) at p = fo, we find 

$' k2 =a 3 a 1 $' k2 Q(k + ), 
where $' k := <J> /+ (/c2 + iO, ik^). The corotating analog of (4.8) is 

0"3$'fc 2 0-l- 



fc 2 

The preceding two equations give 

- vi$'k 2 °iQ = o-i& k2 Qai. 



Premultiplying both sides by — <&' k a\, we find (5.18). 



(5.19) 

□ 



Propositions 5.1 and 5.2 imply the following result. 



Proposition 5.3 Suppose that f' = on the disk. Let B, B, and A(/c) be defined by 



(5.14). Then the spectral functions F(k) and G(k) satisfy the relation 

{B^k^o-iazkB^A+o-xA- 1 ) = -(A+^Al^tB^A^o-io-sAB), keT, 



where A + (k) is defined in terms of F(k) and G(k) in Proposition 4-1 



5.4 The Neugebauer-Meinel disk solutions 

For the BVP denoted by (A) in the introduction, the assumptions of equatorial sym- 
metry and of f' p = on the disk are both valid. In this case the spectral functions F 
and G can be constructed in terms of the known boundary values alone, so that the 
BVP can be effectively solved. The resulting solutions are the celebrated Neugebauer- 
Meinel solutions describing rigidly rotating disks of dust cf. [21. 



Combining Propositions 4.3 and |5.3| we find a 2 x 2-matrix RH problem for the 
function A(k)o-\A~ l (k) in the complex /c-plane with jump across V. In order to 
express our result in the form presented in (21], we formulate this RH problem in 
terms of the 2 x 2-matrix valued function M{k) defined by 



M(k) = a 3 a 1 A(k)a 1 A 1 (k)aia 3 




G 2 (fc)-i v 

F(k) 

-G(k) _ 



(5.20) 



Proposition 5.4 Suppose f is a solution of the BVP denoted by (A) in the intro- 
duction. Let /o := /(-HO) = e 2Uo + ib^ denote the value of f at the origin. Then the 
spectral functions F(k) and G(k) are given by 

F(k) = -M 2 i(k), G{k)=Mn{k), keC, 

where M is the unique solution of the following RH problem: 
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A4(k) is analytic for k € C\T, T = [— ipo,ipo\. 
Across r, M(k) satisfies the jump condition 

S{k)M-{k) = -M + {k)S{k), keT, 



(5.21) 



where M + and M denote the values of M to the right and left of F, respec- 
tively, and S(k) is defined by 



S(k) 



/o/o - 40 2 /c 2 ib + 2iQ,k 
ibo — 2iQk — 1 



k 6 r. 



M has the asymptotic behavior 

M{k) = -a 1 + 0(l/k), k^oo. 



Proof. Defining S by 



S = e 2Uo a Z (Ji{B- v K- l (j Z (TiAB)a z , 



(5.22) 



(5.23) 



(5.24) 



we deduce from Propositions 4.3 and 5.3 that the function Ai defined in (5.20 ) satisfies 



the jump condition (5.21). Evaluating (5.10) at z = and using (5.4), we find the 



following expression for A(fc): 

A(fc) = / + ikQe- 2U °{ai - I)a 3 . 



Substituting this expression for A together with the expression (5.14) for B into 



(5.24), we find that S can be written as in (5.22). The asymptotic behavior (5.23) 



follows from the properties (3.11) of F and G. 



□ 



Remark 5.5 1. The RH problem in Proposition 5.4 coincides exactly with the RH 
problem (2.77) in |21j. This RH problem can be reduced to a scalar RH problem. 



Indeed, the jump condition (5.21) is of the form 

S X M~ = M + S 2 , keT, (5.25) 
where the invertible matrices Si and 52 are simultaneously diagonalizable, i.e., 

S^VDjV- 1 , J = 1,2, 



where T>\ and T>2 are diagonal matrices and V is an invertible matrix. Hence, (5.25) 
can be written as 

v x {v- l Mvy = {v- 1 mv) + v 2 . 

The (11) entry of this equation provides a scalar RH problem which can be solved 
explicitly. In fact, Neugebauer and Meinel were able to combine the auxiliary RH 



problem of Proposition 5.4 with the main RH problem of Proposition 3.4 in such a 
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way that the diagonalization and solution of the combined RH problem yields the 
Ernst potential / directly. 

2. For the BVP denoted by (A), the constant v alue /q = e 2U o-\-ib' of the corotating 
potential /' on the disk is prescribed. In view of (5.8), we find Re(/o) = e 2U ° = e 2U °. 
On the other hand, the value of the imaginary part of /o, which remains unknown in 
the formulation of the RH problem (5.21), disappears in the diagonalization process 

HQ. 



3. The parameter b' can be set to zero without loss of generality, since b' is 
defined only up to an arbitrary integration constant. Hence the solutions of (A) 
are parametrized by the three real parameters Uq, O, and pq. However, we noted 
in Remark |3.3| that F and G do not jump at the endpoints of the contour V for a 
nonsingular solution. Thus, = Ai~ at the endpoints of T and we infer from 
(5.21) that Tr(5(±ipo)) = 0. This imposes one real condition on the parameters 
Uo, Q, and pq: 

|/o| 2 + 4fi 2 pg = l. 

For physically relevant solutions the Ernst potential / should be nonsingular in all 
of spacetime. This imposes further restrictions on the parameters Uq,£1, and po, see 

nana eb. 



A Lax pair singularities 



At the two branch points k = —iz and k = iz of the Riemann surface S z , A = oo 
and A = 0, respectively; thus the Lax pairs (2.2) and (3.1) have singularities at these 
points. In this appendix we analyze in detail the behavior of the eigenfunction <f>(z, k) 
of the linear problem near the branch points. In particular, we find that (p{z, k) is 
nonsingular near these points. Similar statements apply to the eigenfunction &(z, k) 
of the nonlinear problem. 

Proposition A.l For a fixed z = p+iQ, p > 0, the behavior of the map k \— > <f)(z, k) : 
S z — > C near k = —iz, i.e. near the singularity A = oo, is given by 



k — > -iz, (A.l) 
[z, k) near the second 



4>{z, k) = <j){z, -iz) - 2iU z (z)X(z, k){k + iz) + 0(k + iz), 

where <j){z, —iz) is a finite number. The behavior of k i— > 
branch point k = iz follows from (A.l) and the symmetry (2.16). In particular, 
4>{z, ■) is analytic as a map S z — > C near the branch points. 

For a fixed k, Im k ^ 0, the behavior of the map z \— > <f>(z, k) near the branch point 
z = ik is given by 

(p(z,k) = cp(ik,k) + 2U z (ik)X(z,k)(z-ik) + 0(\z-ik\ 3/2 ), z -> ik. (A.2) 



Proof. Fix i £ C with Imfc ^ 0. The Lax pair equations (2.2) imply that there 
exist constants e > and C > such that 



(f) z {z,k) 



k-k 



k+iz U z (z) + hi(z,k), 
$z(z,k) = h 2 (z,k), 



\z — ik\ < e, 



(A.3) 
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where the functions h\(z,k) and li2(z,k) satisfy 

\hi(z,k)\ < Cy/\k + iz\, \h 2 {z,k)\ < C^\k + iz\, \z - ik\ < e. (A.4) 



The equations (A. 3) show that the value (f>(ik, k) of (f)(z, k) at the branch point z = ik, 
if finite, is given by 



(ik, k) = <f>(z, k) + 



,ik] 



k-k 



k + iz 



-U z (z) + h\(z ', k) dz + h,2(z' , k)dz 



Letting z' = ik + pe ld with < p < r 
right-hand side of (A. 5) as the following sum of three terms: 

lk-k 



(A.5) 
\z — ik\ < e. 

ik\, we can write the integral on the 



U z (ik + pe w )e w dp- I h l {ik + pe l ° ' ,k)e Lb 'dp- I h 2 (ik + pe w ,k)e~ w dp. 



o 



ipe 







.10 



-10, 







In view of (A.4), the integral involving h\ satisfies 



rr 

/ h l {ik + pe w ,k)e ie dp 
Jo 



< 



C^fpdp = 2 ^-r 



3/2 



< r < e. 



A similar estimate holds for the integral involving On the other hand, the integral 
involving U z satisfies 



^—^U z {ik + pe w )e ie dp =U z {ik)e ie 
ipe w 



ipe 



dp + 0(r 



3/2^ 



2iU z (ik)^i(k - k)re ie + 0(r 3/2 ), 



0. 



It follows that the integral in (A.5) converges, the value (f)(ik,k) at the branch point 
is finite, and 



(z, k) = <j)(ik, k) - 2iU z {ik)Ji{k -k)(z- ik) + 0(\z - ik\ 3/2 ), 



ik. (A.6) 



Equation (A. 2) is simply an alternative way of writing this expansion which makes 
the choice of bran ches of th e square roots more evident (this choice can be fixed by 
substituting (A. 2) into (2.2) and using that X z = 2 (fc-H Z ) ) • 

We now prove (A.l). For a given z = p + i(, p > 0, the map 4>(z, •) from S z to 
C is analytic in a punctured neighborhood of k = —iz with a possible singularity at 
k = —iz. But equation (A.6) implies that <j)(z, k) is bounded for all z, k near a branch 
point. Since the function cp(z, ■) is bounded near k = —iz, it is in fact analytic in a 
neighborhood of k = —iz. For each z, we may therefore expand 4>(z, •) in a power 
series as follows: 



(z, k) = 4>o(z) + 4>i(z)y/k + iz + 4>2{z)(k + iz) + 



A- 



-iz. 
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Substituting this expansion into the first of the Lax pair equations in (2.2 ), the terms 
of 0(l/y/k + iz) yield 

This leads to the expansion 

<fi(z, k) = cj)(z, —iz) — 2i\J — 2ipU z (z)Vk + iz + 0(k + iz), k — > — iz. 

Equation ( |A,1 ) is an alternative way of writing this equation which makes the choice 
of branches more evident. 



Remark A. 2 The integral representation (2.20) for <f)(z,k) is consistent with the 
expansions ( A.l ) and (A. 2 ). This can be checked by using the following representation 
for U z obtained from (2.21), 



U z {z) 



D{k) 

8ir Jy X(z, k){k + iz) 



rdk. 



B Abel transforms 

In this appendix we use Abel transforms to verify explicitly that the integral represen- 
tation of the solution U (p, £) of the axisymmetric Laplace equation given in Theorem 



2.5 indeed yields the correct boundary values. 
Let U(p,() be defined by 



u(p,0 



i 

4iri 



D{k) 



-.dk, 



(B.l) 



■'Po V(k-() 2 + p 2 ' 

where the branch with positive real part is chosen for the square root, and D{k) is 
given by (2.38) with U(p + i0) replaced by Uo(p), i.e. 

" m U 0p (p) 



( r 1 ' 

£>(*) = -4 I U (0) + \k 2 \ 



V k 2 ~ P' 



;dp 



k = ik 2 , \k 2 \ < po- (B.2) 



We need to show that U = (7o on the disk. 
For z = p + iO just above the disk, we find 



l(=+o" 



'- W k 2 ~ P 2 > 



p <k 2 < po, 
\k 2 \ < p, 
-po <k 2 < -p. 



(B.3) 



Equations (B.l) and (B.3) imply that 

rp 



U(p,+0) 



1 

2^ 



dk 9 



V^fcf 



D(ik 2 ), 0<p<p . 



(B.4) 
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Defining the Abel transform h{k 2 ) of a function h(p) by 



h(k 2 ) 



1 / h(0) , f k2 h p (p)dp 



+ 



7T I k 2 Jo s /]^ - pi 



equation (B.2) can be written as 

D{ik 2 ) = -47r|A&|l7o(|A&|). 



(B.5) 



(B.6) 



Substituting (B.6) into (B.4) and using that the inverse of (B.5) is given by 

h(k 2 )k 2 dk 2 



Hp) = 2 



we infer that indeed U(p, +0) = Uo(p). 

In order to verify the integral representation in terms of the Neumann boundary 
values, we let D(k) be given by (2.39) with U^(p + iO) replaced by U%(p), i.e. 



D(k) = 4 



p0 pUi(p) 

i*2 1 Vp 2 - k 'i 



dp, k = ik 2 , \k 2 \ < po. 



(B.7) 



We need to show that U^(p, +0) = U\{p). We define the following slight variation of 
the Abel transform (B.5): 



h{k 2 



1 d f p0 h{p)pdp 



nk 2 dk 2 J k2 _ jq ' 



(Bi 



whose inverse is given by 



Hp) = 2 / 

Jo 



f° h{k2)k 2 dk 2 



Equation (B.7) yields 



dh 



-D(ik 2 ) = -4TTk 2 U(;(k 2 ) : p<k 2 <p . 



On the other hand, equation (|B.1|) yields 
U C (p, +0) 



1 d 

"4vrdC 



p0 D{ik 2 )dk 2 
(=o+J-p \Jiik 2 - C) 2 + p 1 



(B.9) 



(B.10) 



(B.ll) 



The function g(p, (, k 2 ) defined by 

g(p, (,k 2 ) = i log (k 2 + + W(ik 2 - C) 2 + p 2 ) 

satisfies 

dg_ = 1 

dk 2 v/^-O' + P 2 ' 
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Thus, integrating by parts in (B.ll) and using that D(±ipo) = 0, we find 



UA P , +0) 



1 d 



p0 d 

9{p,C,k2)-jrD(ik 2 )dk2. 
dk 2 



We can now interchange the orders of differentiation and integration. The identity 



dg . dg i 
dC~ l dh~ V(^2-C) 2 + /^' 



implies that 



1 f Po 



^ J-po V(ik2-C) 2 + P 2 \ (=0+ dk ? 



D{ik 2 )dk 2 . 



Using (B.3) together with the fact that -jj^D{ik 2 ) is an odd function, we infer that 



1 f Po 



2tt J \ v /fc| - p 2 dk 2 



D{ik 2 )dk 2 



(B.12) 



Substituting (B.10) into (B.12) and using (B.9), we find that indeed U^(p,0) = U\(p). 
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